* Make use of survey features in GLES data
* overall weights -> pweight, region -> strata, sampling points

svyset vnvpoint [pweight=w_ipfges_1] , strata(distost)

